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Abstract 

We study the problem, initiated in of finding randomized and quantum complexity of initial- value 
problems. We showed in |^ that a speed-up in both settings over the worst-case deterministic complexity 
is possible. In the present paper we prove, by defining new algorithms, that further improvement in upper 
bounds on the randomized and quantum complexity can be achieved. In the Holder class of right-hand 
side functions with r continuous bounded partial derivatives, with r-th derivative being a Holder function 
with exponent p, the e-complexity is shown to be O ((1/ eY^'^^'^P^^^^^) in the randomized setting, and 
O on a quantum computer (up to logarithmic factors). This is an improvement for 

the general problem over the results from The gap still remaining between upper and lower bounds 
on the complexity is further discussed for a special problem. We consider scalar autonomous problems, 
with the aim of computing the solution at the end point of the interval of integration. For this problem, 
we fill up the gap by establishing (essentially) matching upper and lower complexity bounds. We show 
that the complexity in this case is Q ((l/e)^/(''+''+^/^)) in the randomized setting, and Q ((l/e)^/(''+''+^') 
in the quantum setting (again up to logarithmic factors). Hence, this problem is essentially as hard as 
the integration problem. 
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1 Introduction 



Significant progress has been made in recent years in the field of quantum complexity of numer- 
ical problems. Integration ( [pl]| , followed by [^) was the first problem to be so studied. Other 
problems were next analyzed, such as approximation |Q, |5| and path integration [||]. The only 
paper that has studied the randomized and quantum complexity of initial-value problems for 
ordinary differential equations is This paper showed that we can achieve a nontrivial speed- 
up by going from the worst-case deterministic setting to the randomized or quantum settings. 
The idea in was to use the optimal deterministic algorithm based on integral information 

^], and replace integrals in a suitable way by optimal randomized or quantum approximations 

10| , We recall the results from in Theorem 1. 

. n the present paper, we show that further improvement in upper bounds on the randomized 
and quantum complexity is possible. We first define a new deterministic integral algorithm for 
initial- value problems (Section 3). Although this algorithm is not optimal in the determinis- 
tic worst-case setting, it is better suited for randomization and implementation on a quantum 
computer than the algorithm used in Randomized and quantum algorithms are defined by 
a suitable application of optimal randomized and quantum algorithms for summation of real 
numbers [||, ^ (Section 4). The reduction of the total cost is achieved due to a better bal- 
ance, compared to the algorithms from Q, between the deterministic and random components 
of the cost. 

New upper bounds on the complexity are shown in Theorem 2 in Section 5. In the Holder class 
of right-hand side functions with r continuous bounded partial derivatives, with r-th derivative 
being a Holder function with exponent p, the e-complexity is shown to be (up to logarithmic 
factors) O |^(l/e)-'^/('"+^+^/^)^ in the randomized setting, and O (^{l/s)^^''^'^^'^^^'^^^ on a quantum 
computer. Noticeable improvement in both settings is thus achieved, compared to the bounds 
from Theorem 1. The gap between upper and lower complexity bounds is reduced (but still not 
cancelled). 

In order to further reduce the gap between the bounds, we turn to a special case of the general 
problem. We study in Section 6 the complexity of computing the solution of a scalar autonomous 
problem at one single point. In we only showed (non-optimal) upper bounds on the random- 
ized and quantum complexity of this problem. The question about lower bounds was left open. 
We provide essentially matching upper and lower complexity bounds in Theorem 3. Upper 
bounds are established by using a bisection argument, while lower bounds by reducing the prob- 
lem to the summation of real numbers. Up to logarithmic factors, the complexity turns out 
to be G ^(l/e)^/(^+''"'"^/^)^ in the randomized setting, and B ^(l/e)^/^^"'"''^-^)^ in the quantum 
setting. The gap between upper and lower bounds is thus essentially closed. Up to logarithmic 
factors, the problem considered turns out to be as difficult as the integration problem. 

2 Preliminaries 

We deal with the randomized and quantum solution of a system of ordinary differential equations 
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z'(t) 





with initial conditions 

z'{t) = f{z{t)), te[a,h], z{a)=rj, (1) 

where / : R"^ R*^, the initial vector t] is in R*^, and the solution z maps [a,b] into R'^. We 
assume that /(r/) ^ 0. 

This formulation covers nonautonomous systems z'{t) = f(t,z{t)) with / : R'^^-'^ R'^, which 
can be written in the form (||) by adding one scalar equation: 

1 

f{u{t),z{t)) 

with an additional initial condition u{a) = a. We assume that the right-hand side function 
/ = [/^, . . . , Z'^]^ belongs to the Holder class F^'^. Given an integer r > 0, a number p € (0, 1], 
positive numbers Dq, Di, . . . ,Dr and H, we set 

F^^P = {f:K''^K''\ /gC'-(R'^), |aT(2/)[< a, i = 0,l,...,r, 

\d'f{y)-d'f\z)\<H\\y-zr, y,^€R", J = l,2,...,d}, (2) 

where 5*/"' represents all partial derivatives of order i of the j-th component of /, and || • || 
denotes the maximum norm in R'^. We assume that p = 1 for r = 0, which assures that / is 
a Lipschitz function. 

We formulate the problem and shortly recall basic definitions concerning randomized and quan- 
tum settings. Our aim is to compute a bounded function / on [a, b] that approximates the solution 
z. Letting {xj} be the uniform partition of [a,b], so that Xi = a + ih with h = {b — a)/n, we 
will construct / based on approximations ai{f) to z{xi) for i = 0,1,..., n. We assume that 
available information about the right-hand side / is given by a subroutine that computes values 
of a component of / or its partial derivatives. In the randomized setting, we allow for a random 
selection of points at which the values are computed. On a quantum computer, by subroutine 
calls we mean applications of a quantum query operator for (a component of) /, or evaluations 
of components of / or its partial derivatives on a classical computer. The transformation (p that 
computes / based on available information is called an algorithm. 

To be more specific, let (Q, S, P) be a probability space. Let the mappings a; E i-^ of (/) be 
random variables for each / G F^'P. By an algorithm we mean a tuple 

</'=({«0(-),<(-),---,«n(-)Wn,^), (3) 

where -0 is a mapping that produces a bounded function based on ag (/), iiif)-, • • • , d^if)-, 

l-{t)=i^{a^{f),a^,{f),...,a^{mt), (4) 
for t S [a,b]. The error of (/> at / is defined by 

e'^(,/.,/)= sup \\z{t)-l^{t)\\. (5) 

te[a,b] 

We assume that the mapping w £ — > e'^((/>, /) is a random variable for each / G F^'^. 

In the randomized setting, the error of (j) in the class F^''P is given by the maximal dispersion of 

e'--<^(</),F^'^)= sup (Ee-(0,/)2)i/2, (6) 
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where E is the expectation. (We could consider as well the maximal expected value of e'^{(j), /); 
this would only change the constants in our results.) The cost of an algorithm (p in the random- 
ized setting is measured by a number of subroutine calls needed to compute an approximation. 
For a given e > 0, by the e-complexity of the problem, comp''^'^'^(F''''', e), we mean the minimal 
cost of an algorithm (p taken among all (/> such that e^^'^'^{(l), F^'P) < e. 

On a quantum computer, the output of an algorithm is also a random variable (taking a finite 
number of values). The randomness in the quantum setting results from quantum measurement 
operations |^]. The right-hand side function / can be accessed through applications of a quantum 
query operator on a quantum space (defined through values of components of /). Evaluations 
of components of / or its partial derivatives on a classical computer are also allowed. For 
a detailed discussion of the quantum query operator, and of the effect of quantum measurement, 
the reader is referred to The error of an algorithm </> at / in the quantum setting is again 
given by (^), and the error of 4> in the class F^'P by 

gquant^^^ ^r,p^ J) ^ SUp iuf { q| P{ e'^ ((/),/)>« } < 5 }, (7) 

for a given number 5, where < (5 < 1/2. For e > 0, (0) implies that the bound e^{(j), f)<£ holds 
with probability at least 1 — 6 for each / iff e^"'^"*((/>, F^'^, 6) < e. Hence, 1 — is the (minimal) 
success probability in computing an e-approximation. 

The value of 5 is usually set to 5 = 1/4. The success probability can then be increased to be 
at least 1 — 6 (for arbitrarily small 5) by computing component by component the median of 
clog 1/6 repetitions of the algorithm, where c is a positive number independent of 6, see Q. 
The cost of an algorithm (j) in the quantum setting is measured by the number of quantum 
queries, together with the number of classical evaluations of / or its partial derivatives, needed 
to compute an approximation. For a given e > 0, by the quantum e-complexity of the problem, 
comp'^^^^^{F^'P,e,6), we mean the minimal cost of a quantum algorithm (j) taken among all (p 
such that e^'^'^'^*((/),F^''',(5) < e . 

We now recall upper and lower bounds on the randomized and quantum complexity for prob- 
lem (|l|) obtained in (We write below log for log2, although the base of the logarithm is not 
crucial.) 

Theorem 1 (Q) For problem we have that 

... r + p+3/2 \ 

comp-d(F^-P, s) = O ( ( i) '-"-^^'^--^^ log 1 , (8) 



Moreover, for d >2 



comp^"^'^*(F"''',e,d) = ^ 



r+p+2 

1 \ (r+p+1)^ /II 



log - + log - 1 . (9) 
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1\ r + p+l/2 



(10) 
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and, forO<6< 1/4, 

compi"^°*(F'^''',e,5) > compi"^'^*(F^'^', e, 1/4) = n (^(^^^"^'^ • (H) 

The constants in the O- and Q-notation only depend on the class F'^'P, and are independent ^ 
e and 5. 

In the deterministic worst-case setting, if only the values of / or its partial derivatives can be 
accessed, the complexity of problem (|l|) is Q(e~^^^'^~^P^). Hence, Theorem 1 shows a speed-up 
in both randomized and quantum settings over the deterministic setting for all r and p. Note 
also that there is a gap in the randomized and quantum settings between the upper and lower 
complexity bounds given in Theorem 1. 

In this paper, we show that further improvement in upper bounds on the randomized and 
quantum complexities is possible (Theorem 2). We start in the next section by defining a new 
deterministic algorithm that will be used to design randomized and quantum algorithms in Sec- 
tion 4. 

In the next sections we shall need results on randomized and quantum computation of the mean 
of real numbers, which we now recall. Suppose we wish to compute the value 

1 ' 



S = -Y.x„ (12) 

for — 1 < Xj < 1. The e-complexity of this problem in the randomized setting is defined as 
the minimal number of accesses to xi, . . . ,Xs that is sufficient to find a random approximation 
to S with expected error at most e, E\A^ — S\ < e. It is proportional to 

min{s,(l/e)2} (13) 

due to the result of Mathe, see for a discussion Q. Note that E]^'^ — 5*1 < e implies that 

P{|A^ -S\> 4e} < 1/4. (14) 

On a quantum computer we can do better than this. The probabilistic error criterion (p^) is 
used in the quantum setting, and the cost of an algorithm is measured by a number of quantum 
queries (quantum accesses to xi, . . . , Xg). It is shown in [l[ (upper bound) and ||^ (lower bound) 
that the quantum complexity of computing the mean is proportional to 

min{s, 1/e}. (15) 



3 Deterministic Algorithm 

We define a deterministic integral algorithm for solving (|l|), which will be the subject to ran- 
domization and implementation on a quantum computer in the next section. 
Let m,n > 1. Define {xj} to be n + 1 equidistant partition points of [a, b], so that Xj = a + ih for 
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i = 0, 1, . . . , n, where h = {b — a)/n. Let {z*} define a partition of each interval [xi, XiJ^i] with 
m + 1 equidistant points zj = Xj + j/i for j = 0, 1, . . . , m, with /i = (xj+i — Xi)/m. Let t/q = ??. 
By induction, we define sequences {y*} and {y*} as fohows. For a given y* we set = y^. 
Given y*- , by z^j we denote the solution of the local problem 

z',,{t)=f{z,,{t)), z,,{z})=y}. (16) 

Letting l*j{t) be defined by 1*^(1) = i'^/kl)z*l''\zj){t - z^^ for t £ {z),z\^^, we set y)^-^ = 

k=0 

^j{Zj^i) for j = 0, 1, . . . , m — 1. Finally, we define the function /* in [xj, Xj+i] by l*{t) = l*j{t) 
for t G [z*-, and we compute the approximation to z{xi^i) by 

yhi=y:+ J fm))dt (0<i<n-l). (17) 

The approximation I to the solution z of ([l|) in [a, h] is defined by 

l{t)=l*{t) iorte[xi,x^+i]. (18) 

Compared to the algorithm used in Q, the construction above is based not only on the points 
{xi}, but also on the finer partition given by {zj}. The approximation /* in [xj, Xj+i] is computed 
by successive applications of Taylor's method with step size h. 

In the sequel, we shall need an error bound for /* in The following lemma, stated 

without proof, is a standard result for Taylor's method, showing the dependence of the error on 
the length of the interval of integration. Let z* be the solution of the problem 

z'(t) = f{z{t)), t € [x„x,+i], z{xi) = y*. (19) 

Lemma There exists a constant M depending only on the parameters of the class F^'^ ( and 
independent of i, y^ and n) such that 

sup \\z*{t)-l*{t)\\<MhJf+P, 

for sufficiently small h (Lh < In 2, where L is a Lipschitz constant for f). ' 

The algorithm defined above is not optimal in the deterministic worst-case setting. It follows 
from this Lemma and the results from that its worst-case error in [a, b] in the class F^''^ 
is 0{1 / {n{nmY~^P)) . This is achieved by using Q{nm) evaluations. With the same number of 
evaluations it is however possible to get error 0(1 / {nmY'^'''^^ ) , see [Q]. 

In order to define randomized and quantum algorithms, we express (|l7| ) in an equivalent form. 
Defining 

<(y) = i2hf^'\y))(y-yj)' (20) 
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and 

9rM) = ^ [finMi + ^h)) - wt^inM + ^h))) , n G [0, 1], (21) 
we can write ([l7|) as 



m—l -'i'^ m—l "'^ 



?/*+! = / <,(^*,W)rfi + /i^+^+'E mA^)du. (22) 

i=o i=o ;^ 



Arguments similar to those used in the proof of Lemma in Q yield (after replacing the interval 
[xj, Xi+i] by [zj, hhy h and y*,l*,w* by y*-, l*j,w*j, respectively) that the functions gij are 

in C^^\[0, 1]), and the derivatives of gij of order 0,1, ... ,r are bounded by constants depending 
only on the parameters of the class F^'^. Moreover, 

\\9ff{u)-gl^\u)\\<H\u-u\P, u,ue [0,1], 
where H is a constant depending only on the parameters of F^'P. 

4 Randomized and Quantum Algorithms 

We shall denote approximations obtained in randomized and quantum algorithms by the same 
symbols as we did in the deterministic algorithm, omitting only the asterisk. In particular, 
the approximation to z{xi) is denoted by yi. We start with yo = r]. For a given yi we put 
Uo = Hi, and denote by Zij the solution of (|16|) (with the initial value computed for yi). 
We compute lij in a same way as l*j (with y^ instead of y*), and we set yj+i = lij{z^j j^i). 
Approximations U in [xi,rEj+i] are defined to be equal to Uj in each subinterval [zj, Zj_^_i], and 
the polynomial Wij is constructed in the same way as w*j, with y* replaced by yi. 
The approximation at Xj+i is defined by 



m-l 

y,+i = y,+ 5] I Wij{h,it))dt + mh'+P+^Mf), (23) 



m-l 

E 

j=0 

where ^j(/) is a randomized or quantum approximation 



1 m—l . 

A{f)^-J2 9^A^)du. (24) 
i=o i 

The approximation I in [a, 6] is defined by l{t) = /^(t) for t G 

For comparison, in |§ we had m = 1 and j4j(/) was taken to be optimal randomized or quantum 

1 

approximation to the integral / gio{u) du. 
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Here, we define Ai{f) in a different way. Let Qfj{f) be the mid-point rule approximation to 



/ Qijiu) du based on N points, 



N-l 



k=0 



Consider tlie first-stage approximation (without computing it) 



(25) 



-j^ m— 1 „ -j^ m— 1 2 m— 1 

j=o n j=o ]=o k=o 



(26) 



We define Ai{f) to be the optimal randomized or quantum approximation (computed component 
by component) to the right-hand side mean of mN vectors in (|2^). 

Consider first the quantum setting. Let ei > 0. For i = 0, 1, . . . , n — 1, let Ai{f) be a random 
variable such that 

, m-lAf-l N o 

Mf)-^i: E5.,K) <^i}>4 (27) 

for all / G F^'P. To compute Ai{f) it suffices to use of order minjmA'", 1/ei} quantum queries 
for computing each component of the mean, see (|T5|) . (A number of repetitions dependent on d 
is also needed to keep the success probability at least 3/4 when passing from components to 
the vector norm. This changes the cost by a constant factor only.) To increase the success 
probability, we take the median (computed component by component) of k results Ai{f), where 



A: = e log 



1 



1 - (1 - (5)V" 



0(log n + log 1/6) 



(with absolute constants in the Q- and 0-notation). We get a new approximation, denoted by 
the same symbol Ai{f), such that 



i- m~lN-l N 



l/n 



(28) 



j=0 k=0 



This yields that 



m-l N-l 



Mf) 



mN 



E E 9^3 (^k] 



j=0 k=0 



< ei for i = 0,1, 



1 >1 



(29) 



The cost of computing Ai{f) is O (n(log n + log 1/6) minjmA^, 1/ei}) quantum queries. 

In the randomized setting, we compute each component of the mean using the algorithm with 

expected error at most ei/4, and cost proportional to min{mA^, (1/ei)^}, see (|l3|). Inequality 

(|l^ ) then holds with e := ei/4. We next proceed as in the quantum case above to compute 

Ai{f) such that ( p9[ ) holds. For this, we need O (n (log n + log 1/6) min{mA^, (1/ei)^}) function 

evaluations. 
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The deterministic part of the cost of algorithm (23) consists of computing coefficients of kj and 
Wij for J = 0, 1, . . . , m — 1, for which we need cm evaluations of partial derivatives of / of order 
0,1, ... ,r, where c only depends on r and d. The computation of the integrals of Wij does not 
require new evaluations. Taking into account all indices i and j, we need in total cnm evalua- 
tions of / or its partial derivatives. 

5 Upper Bounds on the Randomized and Quantum Complex- 
ity 

We now prove new upper bounds on the complexity of (|l|). 

Theorem 2 For problem there exist constants Pi and P2 depending only on the parameters 
of the class F'^'P such that for sufficiently small e and 5, 

/I X l/(r+p+l/3) 1 

couiv {F"'", e)<Pi[-\ log - (30) 

and 

comp'i-'^*(F^''', E, 5) < P2 i^-j (log - + log - j . (31) 

Proof We analyze the error of the algorithm defined in the previous section. Let = z{xi) — yi. 
Since 

m-1 ^^+1 

z{x,+i) = z{xi) + f{z{t))dt, (32) 

J 



by subtracting (^31) we get that 



m-1 ^^+1 m-1 ''^'+^ 



( lb J. n lib— _L n 

Ci+i =ei+Y. / (/(z(t)) - fikjit))) dt+Y. ifihjit)) - Wijik.m dt - mh^+P+'Aif). 
i=o i j=o i 

J J 

(33) 

Hence, 



m— 1 -'^^ -. m— 1 j, 

ei+i|| < ||e,|| + 5] j \\f{zit))-fikjitmdt + rnh'-'^P+'-J2j9ijiu)du-Aiif) (34) 



i=o 



for i = 0,l,...,n — 1, where the function gij is defined for yi. 

Let Zi be the solution of (ll9|) with the initial condition z{xi) = yi. Using the well known 
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dependence of the solution on initial conditions and the Lemma above, we get for t € [zj,z^j_^i\ 
that 

ii/(z(t)) - fik.m < \\f{z{t)) - fiziim + wfim) - fihm 

< L\\z{t) - zi{t)\\ + L\\zi{t) - kj{t)\\ 

< Leyip{Lh)\\ei\\ + LMkh'^+P 

for Lh < In 2. Inequality together with (^9|) yield now that the inequalities 



|ei+i|| < llcill (1 + /iLexp(/iL)) + LMK^ K'+p 



m-l / 1 ' 

m E \ ! 9ij{u) du - J2 9ij{uk) 
j=0 \0 k=0 . 



(35) 

hold for i = 0, 1, . . . , n — 1 with probability at least 1 — 5. We now take into account the error of 
the mid-point rule, and solve the resulting difference inequality with cq = 0. With probability 
at least 1 — 6, we get that 



<C{h + l/N + ei)h'^+P, i = 0,l,...,n, 



(36) 



for a constant C depending only on the parameters of the class F'^'P. The total cost of computing 
yQ,yi, . . . ,yn is equal in its deterministic part to cnm evaluations of partial derivatives of /. 
The non-deterministic part includes 

O (ra(log n + log 1/6) minjmA^, 1/ei}) quantum queries in the quantum setting, and 
O (ra(log n + log 1/S) minjmA^, (1/ei)^}) evaluations of / in the randomized setting. 
It follows from (p^) with N >n and ei = 1/n that 



Icill < ChhT+P , 



(37) 



for i = 0, 1, . . . ,n, with probability at least 1 — 5 (and a different constant C). Passing to the 
approximation over [a, 6], we get for t S the inequality 



\\z{t) - l{t)\\ < \\z{t) - Zi{t)\\ + WUt) - k{t)\\ < exp(/iL)||ei|| +MhJf+P 
This yields that with probability at least 1 — 6, the error bound 

sup \\z{t)-l{t)\\ < ChhJ'+P 

te[a,b] 



(38) 



holds, with the constant C depending only on the parameters of the class F^'^p. 
Consider the quantum case. Neglecting for a while the logarithmic factors, we have that error 
0{1 / {n{n'mY^P)) is achieved with cost 0{n'm + v?). It is easy to see that the best choice in this 
case is m = n. With a total number of k quantum queries and deterministic evaluations, we 
then achieve the error bound 



sup \\z{t)-l{t)\\ < Cl ^ 



(39) 



with probability at least 1 — 6. This holds for all / € F^'P, and a constant Ci depending only 
on the parameters of the class F'^'P. Hence, to compute an e-approximation I such that 

sup \\z{t) - l{t)\\ < e 
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with probability at least 1 — 6 for each / G F'^'P, the algorithm uses 



O (logl/e + loglA)(l/e 



,l/(r+p+l/2)\ 



quantum queries and deterministic evaluations (the logarithmic factors are again taken into 
account). This completes the proof of Theorem 2 in the quantum case. 

In the randomized setting, we proceed in a similar way, with N > n'^ and m = v? . With k calls 
of / or its partial derivatives (the logarithmic factors are for a while neglected) , we get the error 
bound 

sup \\z{t) - l{t)\\ < C2 k~^'-+P+^/^'^ . (40) 

tG[a,6] 

This holds with probability at least 1 — 5 and a constant C2 depending, as above, only on F'''^ . 
Denote the left-hand side random variable in (|40| ) by , and the right-hand side by h(k). We 
note that 

E(X'^)2= j {X'^fdP{u)+ J {X'^fdP{uj)<K'^5 + h{kf 

X^>h{k) X^<h{k) 

for all / G F'"''', where X is a positive constant, depending only on the parameters of the class 
F'^'P, such that X^ < K. To see that such a constant exists, note that the random variable Ai{f) 
in ([27I) can be assumed bounded by ||^j(/)|| < 2M, where M is a bound on \\gij\\ (otherwise 
Ai{f) = would be a better approximation). Proceeding from (|35|) to (|36|) with ei = 3M, we 



see from (36) that X'^ is bounded (in the deterministic sense) by C}f~^P. Hence, the constant 
K indeed exists. 

Take now k to be the minimal number such that h{k) < e/2, so that k ^ (l/e)i/(^'+P+i/3) ^ and 
set <5 = min{l/2,3eV(4K2)}. Then E (sup^gf^^f,] ||z(t) - /(t)||)^ < for all / G F'^^p, which 

is achieved with cost O ^log(l/e) • (1/e)"'^/^''"'"''"'"^/'^^^ This proves Theorem 2 in the randomized 
setting. 

The upper bounds obtained in Theorem 2 are better than those from Theorem 1 for all r 
and p. For instance, for r = and /) = 1, if we neglect the logarithmic factors, Theorem 1 
gives the bound 0((l/e)^/^) in the randomized setting, and 0((l/e)^/^) in the quantum setting. 
In Theorem 2 the respective bounds are 0((l/e)^/^) and 0((l/e)^/^). Nevertheless, we see from 
lower bounds in Theorem 1 that the gap still remains between the upper and lower bounds. 

Remark 1 

We comment on the proof of Theorem 2, and show a relation to Theorem 1. Looking at (p^) we 

observe that, before starting randomized or quantum computations, we can separate the main 

1 1 1 

part of / gij{u) du by replacing this integral with / Sij{u) du + J {gij{u) — Sij{u)) du, where Sij is 

00 
an approximation to gij. Using / evaluations of gij (/>!), we can define Sij to have the error of 

order l~^'^~^P\ with the cost of one evaluation of Sij independent of I. We can next use randomized 
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1 

or quantum algorithms to compute / {gij{u) — Sij{u)) du. In this way, we get errors ||ej|| of order 



with cost (up to logarithmic factors) 

nm + nml + nmm{mN, (l/ei)'^}, (41) 

where k = 2 in the randomized setting, and k = 1 on a quantum computer. By selecting optimal 
parameters, we get that the minimal (upper bound on the) error achieved with cost k is equal 
to k~^'^~^P^^^^^ in the randomized setting, and k~(^+P+^/'^) on a quantum computer. Hence, by 
admitting / > 1 and by allowing a selection of Sij we do not arrive at better bounds than those 
given in Theorem 2, in which the functions Sij = have simply been taken. 



The upper bounds from Theorem 1 are a special case of (41), and can be obtained for sufficiently 
large N by setting m = 1, ei = n~^^^'^'^'^'^f'^^^ and / = n'^/('^i'+'^P+'^) in the randomized setting, 
and m = 1, ei = n~^^^^^^~^^^ and / = n^/^'''^P^^^ on a quantum computer. 

6 Scalar Autonomous Problems 

In this section, we study the solution of a scalar autonomous problem. The aim is to compute 
the value of the solution at the end point of the interval of integration. We give essentially tight 
upper and lower bounds on the complexity of this problem. In our previous paper Q, no lower 
bounds for this problem were obtained. Upper bounds were discussed together with the general 
problem, which led to weaker results. 

Note that the complexity of approximating the solution at only one single point may differ from 
that of approximating the solution over the whole interval of integration, which is the subject 
of the proceding part of this paper. In particular, upper bounds for the former problem need 
not be valid for the latter one. 

Consider problem (|l]) with d = 1, and the right-hand side function / belonging to the class 

/ E F-.p = F^.p n {/ : \fiy)\>p, yen}, (42) 

for some p > 0. Our aim is to compute the value z{b) with accuracy e by randomized or quantum 
algorithms. Since 

t-a=J^^ds, 
n 

we equivalently look for the solution y* = z{b) of the nonlinear equation H{y) = 0, where 

y 

H{y) = J -j^^ds-{b-a). (43) 

V 
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(The idea of transforming a scalar autonomous problem into a nonlinear equation was exploited, 
for example, in ||2j to derive a class of nonlinear Runge-Kutta methods.) 
Note that ^ ^ 

YtIv -y\< \H{y) - H{y)\ < -\y - y\ , 

for all y and y. 

Given y, the computation of H{y) reduces to the computation of the integral. Suppose that 
we have at our disposal a randomized or quantum algorithm for computing integrals, which 
computes a random approximation A{y) to H{y) such that 

\H{y)-A{y)\<ei (44) 

for some (small) ei > 0, with probability at least 3/4, for any / and y. We denote the cost of 
this algorithm (dependent on a current setting) by c(ei). 

We now define algorithms for computing an approximation to y* = z{b) with error at most e 
with probability at least 1 — 5, for all / E F'^'^. We shall use the bisection method based on 
the values A{y). To get success probability at least 1 — 6, we shall need inequality (H) to hold 
with probability higher than 3/4. Let i* be the minimal index i for which DQ{b— a) / {p2^^^) < ei, 
i.e., i* + 1 = [log {Do{b - a)/{pei))]. We need (H) to hold with probability at least 1 - Si, 
where (5i = 1 - (1 - (5)i/(**+i). To increase the success probability in computing A(y) from 3/4 to 
1 — Ji, we proceed in a standard way by computing the median of k repetitions of the algorithm, 
where 

k = 0(log l/6i) = 0(log {i* + 1) + log 1/5) = 0(loglog 1/ei + log 1/6) . (45) 

Assume that f{y) > (the case /(y) < is analogous). We start the bisection method from 
the interval [ao,/3o] = [r],rj + DQ{b — a)] containing y*, and we set yi = (ao + /3o)/2. Given [aj,/?^], 
we set yj+i = (a^ + /3j)/2 and select the next interval [aj+i,/5j+i] based on the sign of yl(yj+i). 
We stop the iteration at first index i, call it i^^^, for which |j4(?/j)| < 2ei (we shall discuss this 
termination criterion and the correctness of the selection of successive intervals in a while). 
Note that for any j < i* , inequalities 

\H{y,+i)-Aiyi+i)\<ei i = 0,l,...,j (46) 

hold (simultaneously) with probability at least (1 — 5i)* = 1 — 6. 

Assume that ( ^6| ) is satisfied. We show that the number of bisection steps satisfies i"^'^ < i* + 1. 
Suppose that the termination condition is not fulfilled by the i*-th step, i.e., |j4(yj)| > 2ei for 
i = 1,2, . . . ,i* . Then the selection of the interval [ai*,/3i*] made on the basis of A{yi*), as well as 
the selection of all proceeding intervals, is correct. (In fact, it suffices for this that |^(yi)[ > ei, 
since the signs of A{yi) and H{yi) are then the same by (|4^). ) Hence, we have 

- y*\ < \a,* - A.|/2 = \ao - /3o|/2'*+\ 

and 

\Aiyi*+i)\ < \H{yi,+i)\ + \A{yi,+i)-H{yi,+i)\ < |y,.+i-y*|/p+ei < \ao-(3o\/{p2'*+^)+ei < 2ei. 
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Since the termination condition is now satisfied, we have in this case that i = i* + 1. In any 
case, the desired bound on i^'^ holds, as claimed. In terms of £i, we have that 

i'^''<\log{Doib-a)/{peim . (47) 

Take now ei = e/(3Do)- Then, terminating after i^^^ steps, we arrive at the e- approximation 
T/jbis to y*, since 

ly^bi^ - y*\ < Do\H{yi^^.)\ < Do{\A{y,^^s)\ + \H{y,^^s) - A(y,M.)|) < 3Doei = s . (48) 
As in the case of (|46|), this holds with probability at least 1 — 6. 

Summarizing, the described algorithm returns the approximation (random variable) yjbis such 
that the bound |yjbis — z(b)\ < e holds with probability at least 1 — S, with total cost 

c{ei)k i^'^ = 0{c(ei) (log log 1 /e + log 1 /6) log 1 /e) . (49) 

Using known results on integration, we now estimate c(ei). Since / G F^''P, the function l//(s) 
is in the Holder class F^'^ (over the finite interval s € [ri,ri + DQ{b — a)]) with certain parameters 
Dq, Di, . . . , Dr, H depending on Dq, Di, Dr, H and p. 

Consider the quantum setting. There exists an algorithm for computing integrals of l//(s) 
with cost c(ei) = O quantum queries, see This leads to the following 

complexity bound for our problem 

comp'5^^°*(F"'^e,5) = O ((l/e)^/("+''+^) (log log 1/e + log 1/6) log 1/e) . (50) 

Remark 2 

To establish the cost of an algorithm in the quantum setting, we have to count the number of 
applications of a quantum query operator Qf for /. Calculating the cost c(ei) above we have 
taken into account the number of queries Qi/f for 1//. However, a query for 1// for / G F^'^ 
can be simulated by a query for / (and vice versa), see Lemma 4 in [^]. Hence, the upper bound 
in terms of both units remains the same. 

Consider the randomized setting. There exists an algorithm approximating integrals with 
the mean square error (^ bounded by ei/2, and cost c(ei) = O (^{l/eiy^^'^^P^^^'^^^ evalua- 
tions of /. We use it to compute an approximation A{y) to H(y), for a given y. By the Markov 
inequality, error bound ( ^4|) holds for A{y) with probability at least 3/4. We now follow the steps 
between relations (^) and ( ^9| ) above to get the approximation y^bis to z{b) such that 

X":=|y,bis-z(6)| <e (51) 

with probability at least 1 — 6, for all /. By (^9|), the cost of computing y^bis is 

O ((l/e)i/(''+''+i/2) (log log 1/e + log 1/6) log 1/e) . 

To estimate the mean square error of i/jbis, we proceed in a similar way as we did in the final 
part of the proof of Theorem 2. We replace e in ( |5l| ) by e/2, which influences the cost only by 
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a constant factor, and we write 

E(X'^)2= J {X'^fdP{uj)+ J {X'^fdP{Lu)<K^6 + e'^/A 

X">e/2 X'^<e/2 

for all / G F^'^. Here, X is a positive constant depending only on the parameters of the class 
F'''P such that < K. The choice 6 = 3e'^/{AK'^) gives the bound 

sup (E(X'^)2)l/2 < £ ^ 

which is achieved with cost 

O (iQgiog i/e + iog l/e)log 1/e) = O (j^g i/^)2j ^ 

This yields an upper bound 

^Qj^prand(^r-,p^ e) = O ((l/s) ^('^+''+1/2) (log 1/e (52) 

on the complexity. Hence, up to logarithmic factors, we are able to solve our problem at cost of 
one single integration. 

We now turn to lower bounds on the randomized and quantum complexity. Let (p be any algo- 
rithm based on evaluations of / or its derivatives at possibly random points in the randomized 
setting, and on quantum queries for / and classical evaluations of / or its derivatives in the quan- 
tum setting. Assume that <j) computes an approximation to z{b) with error at most e, for any 
scalar problem (|l|) with / £ F^'^. We estimate from below the number of evaluations (queries) 
used by 4>, by reducing the problem to the summation of real numbers. 

Without loss of generality, let [a, b] = [0, 1]. For n > 1, let Aq, Ai, . . . , A„_i be numbers of at 
most unit absolute value, and define the function g{y) = 1/ f{y) as follows. 

Consider the uniform partition of [r], rj + 1/2] with points Vi = r} + i/[2n) for i = 0, 1, . . . , n. We 
let hi G F^'P, where i = 0, 1, . . . , n — 1, be functions with the following properties: 

hi has support , hi{y) > 0, 

max/ii(y) = hi{{yi + yi+i)/2) = ci(yi+i - yif^'', 

Vi+i 

hi{y) dy = C2{yi+i - yiY^''^^, 

Vi 

where ci and C2 are known positive constants depending only on the parameters of the class 
F'^'P (and not on i and n). Such functions are often used in proving lower bounds and their 
construction is well known. 

n-l 

We define g{y) = 1 + ^ihi{y)- Then g £ F'^'^ for sufficiently large n, and the same holds 

1=0 
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for / (with different constants). Since 3/4 < f{y) < 3/2 for sufficiently large n, we liave tliat 
ry + 3/4 < z(l) < r/ + 3/2, and we can write 



z{l) v+1/2 z{l) 



n-1 



J l/f{y)dy= J l/f{y)dy+ J l//(y) = csn-^'^+^+i) ^ A, + z(l) - 77 • (53) 
This yields that 



1 l-z{l)+rj 
where C3 is a known positive constant. 

Consider first the quantum setting. Let C (a random variable) be an e-approximation to z{l) 
computed by the algorithm cp for the right-hand side / defined above. We have that 

\z{l)-C\<e 

with probability at least 3/4 (we take 6 < 1/4). Hence, 



Ci = . , ^ is an approximation to S = — A 



with error at most £2 = £n^~^P/c^ and probability at least 3/4. The lower bound of Nayak 
and Wu, see ([l5|), gives that the number of queries for A = [Aq, Ai, . . . , A^i-i] must be at least 
r2(min{n, 1/62}). This is also a lower bound on the number of queries for 1// (and for /) needed 
in the algorithm 0. We now take n x (l/e)^/^^"'"''"'"^), and we conclude that 

comp^^^°*(F^'^,e,l/4) = ((l/e)i/(^'+''+i)) . 

In the randomized setting, let (E|z(l) - C?)^'"^ < Then (ElS* - CiH^''^ < en''+''/c3 = £2 
and the same inequality holds for ElS" — Ci|- Due to (p^), the number of accesses to Aj must 
be at least Q (min{n, (1/^2 )^})- This is also a lower bound on the number of the number of 
evaluations of / or its derivatives. We now take n x (l/e)i/(''"+P+i/2) ^ get 

comp'^^°'^(F'''^,e) = n ((l/e)i/(^+''+^/2)^ _ 
We have shown the following 

Theorem 3 Consider the scalar autonomous problem described in the beginning of this section, 
with a right-hand side f in the class F'^'P . There exist positive constants Pi (i = 3,4, 5, 6j 
depending only on the parameters of the class F^'P such that, for sufficiently small e and 6, 
the following complexity bounds hold true. 
In the randomized setting 

, . / l/(r+p+l/2) . ^.2\ 

comp'-^\F''P,e)<P3U-] Mog-j (55) 
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and 

, . / -, N l/{r+p+l/2) 

comp"''"^{F'^'P,e)>PA-) . (56) 



In the quantum setting 

comp^"-* {F^'P, e, S) < P, U-j (log log - + log - j log - j (57) 

and, forO<6< 1/4, 

^Q^pquant^^r.p^^^^) > COmp<5^^°* (F^'^, £, 1/4) > Pg f - J J • (58) 



Note that in both randomized and quantum settings upper and lower bounds in Theorem 3 are 
matching, up to logarithmic factors. The question of finding matching upper and lower bounds 
for the general problem still remains open. 
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